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Abstract  We  recently  introduced  SPOC  (Spatial  Processing:  Optimized  and 
Constrained)  [1]  as  an  innovative  replacement  for  the  spatial  spectrum  estimation 
traditionally  done  by  beamformers.  This  work  extends  the  method  described  in 
[1]  by  improving  the  formalism  and  presenting  results  with  at-sea  sonar  data. 

Beamforming  has  traditionally  been  performed  by  conventional  (CBF)  or  adaptive 
(ABF)  processing.  These  both  provide  estimates  of  spatial  spectrum  in  a  given 
time  interval  (snapshot).  A  trade-off  between  detection  and  resolution  governs 
CBF  performance.  ABF  offers  improvement  by  lessening  the  effects  of  interfering 
signals  at  each  beam.  This  is  done  with  an  estimated  covariance  matrix  which 
involves  a  trade-off  in  selecting  the  averaging  interval.  Long  averaging  times 
produce  better  estimates  of  the  true  covariance,  but  shorter  intervals  decrease 
the  harmful  impact  of  target  dynamics. 

SPOC  produces  spatial  spectrum  estimates  without  using  covariances.  It  considers 
an  array  of  N  sensors  and  a  finite  but  large  number  of  sources  of  energy  L 
(L»N)  spanning  the  entire  Direction  of  Arrival  (DOA)  space  (-1<cos0<1).  SPOC 
is  a  constrained  optimization  approach  to  estimating  the  signals  at  each  of  the 
L  DOAs  (or  "beams").  The  complex-valued  signals  are  constrained  so  that  the 
superposition  of  the  weighted  steering  vectors  at  each  DOA  plus  the  estimated 
noise  equals  the  measured  data.  Many  sets  of  signals  generally  satisfy  the 
constraint  so  the  minimum-energy  set  of  signals  is  chosen.  Mathematical  details 
will  appear  in  the  fuil-length  paper. 

As  a  preliminary  demonstration  of  performance  improvement  over  ABF,  we 
compared  APB-98  passive  sonar  data  (system  ABF)  versus  the  same  data 
processed  using  SPOC.  The  processing  was  done  on  every  frequency  bin  at 
the  system  resolution  and  broadband  processing  done  on  both  using  the  same 
standard  "SPED  CS"  algorithm  [2].  The  passive  broadband  grams  demonstrate 
that  more  contacts  are  visible  with  SPOC  processing  than  with  ABF. 

[1]  R.  Bethel,  B.  Shapo,  and  H.  Van  Trees,  "Single  Snapshot  Spatial  Processing:  Optimized 
and  Constrained,"  Sensor  Array  and  Multichannel  Signal  Processing  Workshop  2002  (SAM  2002) 
Proceedings,  pp.  508-512. 


[2]  M.  Bono,  R.  Bethel,  P.  McCarty,  and  B.  Shapo,  "Subband  Energy  Detection  Methods  in  Passive 
Array  Processing,"  MIT  Lincoln  Laboratory  ASAP  Workshop  Proceedings,  2001. 
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ABSTRACT 

The  purpose  of  a  radar  or  sonar  system  is  to  detect  targets  of 
interest.  Beamforming  plays  a  pivotal  role  in  this  task.  Spatial 
Processing:  Optimized  and  Constrained  (SPOC)  is  an 
innovative  front-end  beamforming  method  that  operates  on 
array  sensor  data  as  an  alternative  to  conventional  or  adaptive 
beamforming.  It  uses  a  maximum  a  posteriori  approach  to 
develop  a  constrained  optimization  solution.  A  parametric 
representation  provides  the  framework  for  the  solution 
presented  here.  We  develop  a  recursive  solution  using  past 
history.  This  method  provides  significantly  improved 
resolution  over  conventional  and  adaptive  beamforming. 
Three  real  data  sets  demonstrate  these  performance 
improvements. 


1.  INTRODUCTION 

We  recently  introduced  SPOC  (“Spatial  Processing: 
Optimized  and  Constrained”)  [1]  as  an  innovative 
replacement  for  the  spatial  spectrum  estimation 
traditionally  done  by  beamformers.  This  work  extends  the 
method  described  in  [1]  by  improving  the  formalism  and 
presenting  results  with  at-sea  sonar  data. 

Beamforming  has  traditionally  been  performed  by 
conventional  (CBF)  or  adaptive  (ABF)  processing.  These 
both  provide  estimates  of  spatial  spectrum  in  a  given  time 
interval  (snapshot).  A  tradeoff  between  detection  and 
resolution  governs  CBF  performance.  ABF  offers 
improvement  by  lessening  the  effects  of  interfering  signals 
at  each  beam.  This  is  done  with  an  estimated  covariance 
matrix,  which  involves  a  tradeoff  in  selecting  the  averaging 
interval.  Long  averaging  times  produce  better  estimates  of 
the  true  covariance  but  shorter  intervals  decrease  the 
harmful  impact  of  target  dynamics. 

SPOC  produces  spatial  spectrum  estimates  without  using 
covariances.  It  considers  an  array  of  N  sensors  and  a  finite 
but  large  number  of  sources  of  energy  L  ( L»N)  spanning 
the  entire  Direction  of  Arrival  (DOA)  space  (-1<cos(0)<1). 
SPOC  is  a  constrained  optimization  approach  to  estimating 
the  signals  at  each  of  the  L  DOAs  (or  “beams”).  The 
complex-valued  signals  are  constrained  so  that  the 
superposition  of  the  weighted  steering  vectors  at  each 
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DOA  plus  the  estimated  noise  equals  the  measured  data. 
Many  sets  of  signals  generally  satisfy  the  constraint  so  the 
minimum  energy  set  of  signals  is  chosen. 

As  a  demonstration  of  performance  improvement  over 
ABF  and  CBF,  we  compared  passive  sonar  data  versus  the 
same  data  processed  using  SPOC.  The  processing  was 
done  on  every  frequency  bin  and  broadband  processing 
done  on  both  using  the  same  standard  algorithm  [2],  The 
passive  broadband  displays  in  the  results  section 
demonstrate  that  more  contacts  are  visible  with  SPOC 
processing  than  with  ABF  or  CBF. 


2.  APPROACH 

We  consider  the  sensor  array  problem  with  complex 
narrowband  data.  We  assume  that  the  spatial  energy 
impinging  upon  the  sensor  array  emanates  from  a  dense 
grid  of  far  field  independent  sources,  each  radiating  a 
narrowband  signal.  In  the  limit,  the  set  of  discrete  sources 
approaches  a  continuum.  This  is  a  common  representation 
of  spatial  energy  [2],  The  goal  is  to  estimate  the  energy  of 
each  discrete  source. 

For  implementation  purposes,  we  consider  a  uniformly 
spaced  linear  array  of  N  sensors.  We  assume  a  finite  but 
large  number  of  sources  of  energy  L  ( L»N )  uniformly 
spanning  the  entire  Direction  of  Arrival  (DOA)  space 
u  (-1  <  u  <  +1)  ■  Each  sensor  n  (n  =  1, ...  JV)  receives  the 
sum  of  all  narrowband  signals  radiated  by  sources  /  at 
DOA’s  u/  (I  =  1  for  all  L  sources.  This  observation 
model  is  expressed  as: 

x  =  (1) 

i=i 

where  x  is  the  /Vx  1  complex  narrowband  sensor  data,  y(ui) 
is  the  /Vx  1  array  manifold  vector  for  DOA  uj ,  V(u)  is  the 
NxL  array  manifold  matrix,  f  is  the  scalar  complex-valued 
/th  signal,  and/ is  the  Lxl  signal  vector. 

\fi  I2  is  the  signal  power  at  the  /th  discrete  DOA.  To 
accomplish  our  goal,  we  need  an  estimate  of  these 


quantities.  This  first  requires  estimates  of  /.  With  no 
knowledge  of  observed  sensor  data  x,  the  signal  vector  /  is 
expressed  in  the  a  priori  Probability  Density  Function 
(PDF)  p(  /  ),  which  is  assumed  zero-mean  complex 
Gaussian: 


where  (»)H  represents  the  conjugate  transpose  operation, 
and  Qf  is  the  LxL  signal  covariance  matrix  E[//H  ].  Under 
the  assumption  of  independent  signals,  C  /  =  diaglO/2), 
where  O/2  is  the  expected  power  of  the  /th  signal,  which  is 
unknown. 

We  use  observed  sensor  data  x  to  improve  knowledge  of 
signal  vector  £  via  the  a  posteriori  PDF  p(  £  I  x  ).  In 
particular,  we  find  the  maximum  a  posteriori  (MAP) 
estimate  of  signal  vector  £  by  maximizing  p(  £  I  x  )  with 
respect  to  /with  .y  known.  The  MAP  signal  power  estimate 
of  each  \£  I2  is  computed  from  the  MAP  estimate  of  fi  . 
This  is  accomplished  by  Bayes’  theorem: 

P(/  I  x)  =  Ap(x  I  /)  p(/) .  (3) 

Scalar  A  is  a  normalizing  constant  and  does  not  affect  the 
computation  of  MAP  estimates.  The  conditional  PDF  p(xl 
/)  is  required  in  (3).  We  note  that  (1)  must  be  satisfied. 
Therefore,  we  must  restrict  signal  vector /to  values  which 
meet  this  requirement.  The  implication  in  (3)  is  that  the 
conditional  PDF  p(xl/)  is  zero  for  any  signal  vector  / 
which  does  not  satisfy  the  observation  model  (1),  and 
conditional  PDF  p(xl/)  is  equally  likely  for  any  signal 
vector/which  does  satisfy  (1). 

Therefore,  to  find  the  MAP  estimate  of  signal  vector  /  we 
must  maximize  the  a  priori  PDF  p(/)  in  (2)  subject  to  the 
constraint  that  (1)  is  satisfied. 

Substituting  this  result  into  (4)  yields: 

Cost  Function  :  T  =  fHQ^f  (6a) 

Constraint :  x  =  V  f  ,  (6b) 

where  Qf  is  the  theoretical  diagonal  covariance  matrix  for 
the  signal  vector  /  and  the  constant  term  in  (4)  has  been 
neglected.  Also,  the  expression  (4)  has  been  negated, 
resulting  in  a  cost-minimization  problem. 

Then,  the  complex  signal  vector/ which  minimizes  the  cost 
r  in  (6a),  subject  to  the  constraint  in  (6b),  is  the  MAP 
estimate.  Signal  powers  I/I2  for  all  DOA’s  it;  are  the 


desired  output  and  are  the  MAP  estimates  for  signal 
powers. 

The  spatial  spectrum  is  defined  as  the  Fourier  transform  of 
the  stationary  spatial  autocorrelation  function.  We 
emphasize  that  SPOC  is  not  a  spatial  spectrum  estimator  in 
this  sense.  SPOC  also  differs  greatly  from  CBF  and  ABF 
in  that  it  is  not  a  spatial  filter.  It  does  not  require  an 
estimated  covariance  matrix  as  does  ABF.  Also,  SPOC  is 
neither  a  source  detection  technique  [4],  nor  a  DOA 
estimation  algorithm  [5,  6,  7],  The  results  in  Section  4 
demonstrate  the  performance  characteristics  of  SPOC. 

3.  IMPLEMENTATION 

The  desired  quantity  in  (6b)  is  the  vector  of  signals  /.  The 
solution  to  (6a)  and  (6b)  is: 

f  =  CfVH(VCfVH+a2I)-lx  .  (7) 

where  CT  is  diagonal  loading  for  numerical  stability.  The 
desired  output  is  the  signal  vector  power  y  ,  on  a  per- 

scan  basis.  Qf  ■  An  estimate  of  the  diagonal  matrix  Cy  is 
required  in  (7).  This  is  obtained  by  the  computing  the 

diagonal  of  as  an  exponentially  averaged  version  of  f  . 

To  take  dynamic  conditions  into  account,  a  smoothing  operation 
can  be  applied  over  DOA.  In  practice,  we  have  applied  a  simple 

l  ~l2 

FIR  filter  over  beams  to  the  signal  vector  power  /  .  This  is 

roughly  analogous  to  the  assumption  that  the  changes  in  the 
signal  DOA  values  follow  a  first-order  Markov  process. 


4.  RESULTS 

We  implemented  this  algorithm  on  developmental  testbed 
with  real  passive  acoustic  data.  The  results  are  fair 
comparisons  between  SPOC  and  CBF  or  SPOC  and  a 
beamspace  MVDR  ABF.  The  single-frequency  SPOC 
processing  described  above  is  implemented  independently 
for  each  frequency  bin  over  a  broadband  bandwidth.  A 
standard  broadband  processing  technique  results  in  the 
bearing  versus  time  data  displays  below.  Energy  is 
intensity  modulated  where  white  indicates  stronger  energy. 
The  horizontal  axis  corresponds  to  bearing  (DOA)  and  the 
vertical  axis  represents  time.  Three  data  sets  are 
presented. 

These  figures  all  demonstrate  that  contacts  missing  in  both 
ABF  and  CBF  are  detected  in  SPOC.  This  is  a  significant 
improvement  in  beamforming  performance. 


DATA  SET  1 


DATA  SET  2 


DATA  SET  3 


CBF 


SPOC 


Abstract 


Back  to  TOC 


Next  Abstract 


5.  REFERENCES 

[1]  Roy  Bethel,  Benjamin  Shapo,  and  Harry  Van  Trees, 
“Single  Snapshot  Spatial  Processing:  Optimized  and 
Constrained,”  Proceedings  of  the  2002  SAM 
Workshop,  2002. 

[2]  Harry  L.  Van  Trees,  Optimum  Array  Processing: 
Detection,  Estimation,  and  Modulation  Theory,  Part 
IV.  New  York:  Wiley,  2002. 

[3]  J.  Capon,  “High  Resolution  Frequency-Wavenumber 
Spectrum  Analysis,”  Proceedings  of  the  IEEE ,  vol. 
57,  pp.  1408-1418,  Aug.  1969. 

[4]  M.  Wax  and  T.  Kailath,  “Detection  of  Signals  by 
Information  Theoretic  Criteria,”  IEEE  Trans.  Acoust. 
Speech,  Signal  Processing,  vol.  ASSP-33,  No.  2,  pp. 
387-392. 

[5]  R.  O.  Schmidt,  “Multiple  Emitter  Location  and  Signal 
Parameter  Estimation,”  RADC  Spectrum  Estimation 
Workshop,  Griffith  AFB,  pp.  243-258,  1979. 

[6]  R.  Roy,  A.  Paulraj,  and  T.  Kailath,  "ESPRIT  -  A 
Subspace  Rotation  Approach  To  Estimation  Of 
Parameters  Of  Cisoids  In  Noise,"  IEEE  Trans.  Acoust. 
Speech,  Signal  Processing,  vol.  ASSP-34,  No.  4,  pp. 
1340-1342,  Oct.  1986. 

[7]  M.  Wax,  “Detection  and  Localization  of  Multiple 
Sources  via  the  Stochastic  Signals  Model,”  IEEE 
Trans.  Signal  Processing,  vol.  39,  No.  11,  pp.  2450- 
2456,  Nov.  1991. 


